Mean-field Model beyond Boltzmann-Enskog Picture for Dense Gases 
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I propose an extension to Boltzmann BGK equation for Hard Spheres. The present model has 
an -ff-theorem and it allows choice of Prandtl number as an independent parameter. I show that 
similar to Enskog equation this equation can reproduce hydrodynamics of hard spheres in dense 
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^ . INTRODUCTION 
•-^ . 

' An overwhelming majority of fluid flow problems of physical and engineering interest cannot be solved using micro- 
scopic simulation methods, such as molecular dynamics, due to the enormous number of degrees of freedom constituting 



the macroscopic systems. In such a scenario, mesoscale descriptions in terms of one particle distribution function, 
such as Boltzmann equation, provide important tools for understanding transport phenomena beyond phcnomenolog- 
' ical hydrodynamic descriptions in terms of Navier-Stokes-Fourier equations. Indeed, the nonlinear Boltzmann kinetic 
equation can accurately predict a wide range of physical properties and flow profiles for low density gases even in 
I ' states very far from equilibrium (see for example 

However, technical difficulties encountered in solving (analytically or numerically) the Boltzmann equation, a non- 
^ ' linear integro-differential equation for the time dependent distributions in a six dimensional phase space, limits its 
application in practice. During the last few decades, this technical problem has been solved for the Boltzmann equa- 
tion in two very important regimes. Firstly, for highly non-equilibrium situations associated with supersonic flows (in 
general for high Mach number flows), direct simulation Monte Carlo (DSMC) method was applied with remarkable 
success (see for reviews dH). Secondly, for very low Mach number flows lattice Boltzmann method is remarkably suc- 
cessful in both hydrodynamic regime as well as transitional regime (see for example 0, 0, 0] ) ■ The lattice Boltzmann 
q method rely on an approximate form of Boltzmann collision term known as Bhatanager-Gross-Krook (BGK) collision 
, approximation. The model Boltzmann equation with BGK collision term retains almost all qualitative features (such 
l— — '■ as correct conservation laws, fl-theorem) of the Boltzmann equation Indeed, the BGK model can be classified as 
\ the first truely successful phcnomcnological model at the level of one particle distribution. The mathematical sim- 
^ ' plicity of this model is often used to obtain exact and semi-exact analytical solutions which can help to understand 
^\ , the hydrodynamics well beyond Navier-Stokes equations [H, i|. The strength and limitations of this model along with 
ways to make it quantitatively accurate (without destroying the basic features such as ii-theorem) is well understood 

^: UMMM- 

• ■ For hard spheres, in order to describe the fluid transport in dense regimes, Boltzmann equation was extended by 
Enskog and further modified by van Beijeren and Ernst (known as the revised Enskog theory (RET)) [HI, HI]. Similar 
qq ■ to Boltzmann's model of diluete gas, particles motion in these models is decomposed into two parts: propogation at 
\ constant velocity followed by collisions in which exchange of momentum between particles happens. However, unlike 
^ ■ Boltzmann model, collisions are understood to be non-local events due to the presence of finite size particles. This 
idea of non-local collisions due to finite size of particles behind Enskog or RET extension of Boltzmann equation 
is borrowed from Van-der Waals' picture of excluded volume in dense gases. However, mainly due to the non-local 



X 



collisions, Enskog extension of Boltzmann model leads to even more intractable form of nonlinear integro-differential 



equation. Thus, it is not surprising that it took almost fifty years to prove even the existence of H -theorem [I3j and 
so far only modest engineering and physical application of dense fluid is modeled via Enskog equation. The non-local 
nature of collision is difficult to handle both for Monte-Carlo method as well as for kinetic modeling via simplified 
phenomenological theories of BGK type. 

From engineering perspective, a dense gas model, where non-ideality can be added as extra terms over a rarefied 
gas model (either of Monte-Carlo type or BGK type), is an extremely desirable solution. In case of DSMC model 
of Boltzmann equation, Alexander et. al proposed a simple modification of propagation step which gave the correct 



equation of state but failed to reproduce the Enskog transport coefficients |14| . An important progress for such 
modeling approaches was reported in Ref. 15, 3], where BGK-like collision terms for Enskog equation were proposed. 
It needs to be remarked that most multiphase extensions of lattice Boltzmann method uses simplifications pertinent 
to incompressible flow in Dufty et. al model 15] to mimic hard-sphere repulsion (see [f| for details of multiphase 
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lattice Boltzmann method). The main idea behind these works was to compute the effect of Enskog collision term on 
momentum and energy balance and explicitly add it in momentum and energy balance equation as a correction to 
the BGK collision term. This approach gave correct viscosity coefficient and equation of state but there is neither H 
theorem nor correct thermal conductivity for these type of models. 

The goal of the present work is to fill this gap and construct a phenomenological model of fluid transport at 
the meso-scale level. Similar to hydrodynamic description given by Navier-Stokes equations, I demand that a good 
phenomenological theory at mesoscale should fullfil the following criteria 

• It must obey conservation laws and second law of thermodynamics(i?-theorem). 

• It must reproduce correct thermodynamic equation of state. 

• It must reproduce correct Enskog transport coefficients in hydrodynamic regime. 

The basic starting point for such modeling is exact dynamics of one particle distribution function / (defined as 
probability of finding a particle at location x with velocity v at time t) as a function of phase variable z = (x, v), 
given by BBGKY hierarchy as 

$/(*!,*) +0 a [/(Zl,t) Vi] = -Ai/a(B lj Z 2j t) (1) 

where /2(1, 2) denotes the two particle distribution function and A\ is an operator whose exact form is not important 
for the present discussion (For details of this equation see for example [l7j])- an d the hydrodynamic fields {p,j,T} are 
defined in terms of one particle distribution function / as 

A closed form kinetic equation is obtained if I set A\ f2 — »/(/), where J ^ typically collision integral, maps functions 
onto functions. Thus for example, ji = /2(1,2,/) with explicit time dependence of /2 entirely contained in / was 
proposed by Bogoliubov to derive the Boltzmann equation. Traditional mechanistic view-point starting from BBGKY 
hierarchy is that non-ideality coming from repulsive part of molecular forces needs to be modeled via three-particle 
and other higher order collision and free-flight of molecules remains unaffected by it. However, so far only widely used 
model is mean-field approximation of Enskog, where all higher order effects are lumped in non-local collisions. In 
the present work, I am advocating an alternate mean-field picture of non-ideality in dense system. According to this 
picture in dense system, dominant change is alteration of free flight of a hard sphere particle. The key new ingradicnt, 
I propose is that A\ f% will also modify the free propogation step of the Boltzmann type equations. In particular, I 
propose to write generalised kinetic equation for f(z,t) as 

d t f(z,t) + d a [f(z,t)v]=J (3) 

where v is some unknown propagation velocity. The physical picture behind such a modification of free propagation can 
be understood using Bogoliubov hypothesis for the equilibration of a non-equilibrium gas, which assumes separations 
of time scale by considering mean free time r of a molecule to be much larger than mean time spent in the interaction 
domain of the another molecule [13] ■ In dilute region, where reduced density is small, using this hypothesis he 
derived Boltzmann equation. Now for hard-sphere in dense region the dynamics in mean free time r has to be more 
complex due to collective effects. In this time scale, for dense system we also need to consider the effect of collective 
motion (hydrodynamics) on individual particles. Such corrections from local collision can also be incorporated by 
analyzing the short time motion of tagged particle due to entropic force generated by particle-particle correlations and 
hydrodynamics force generated by other particles. In other words, apart from usual hydrodynamic forces the tagged 
particle might also experience effective forces, which have purely entropic origin. Unlike ring kinetic theory approach, 
I am not trying to derive such an equation from the first principle but trying to model it based on physical intuitions 
and limiting behaviors of the system. Firstly, we know the two limiting behaviors of such a system in rarefied and 
extremely dense regime quantified by compressibility factor x defined as 

X=^—l = — O-P— 1 (4) 
X P RT P R\ 9 dp ) ' ( ' 

where the excess entropy s nld (p) and pressure p(p, T) as a function of density p and temperature T is known from 
equilibrium statistical mechanics (For example Van-dcr-Waals or Carnhann-Starling approximation). In a rarefied 
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system (% — > 0), present model should recover Boltzmann description of free propagation step (v = v). while for 
extremely dense system (x — * oo), all hard-spheres should pack together and move with the collective velocity of the 
system, which means v = v = u. 

In the present work, I ask the question what is the most general form of the propagation velocity v (which can be 
a function of the moments of /) in Eq. Q. In order to do so, first we need to recognize that the propagation velocity 
can be written as a formal Hermite expanion v = v + Y^Lo a ^ H*-™^), in terms of dimensionless peculiar velocity 
£,a = ("q — u a ) /v2 RT. Furthermore, we need to recognize that the condition of having correct continuity equation 
(obtained by integrating Eq.([3]) over v) itself severely restricts the choice of v. For example, a( ) = and a^™) = 
for all n > 2. For n — 2 only the trace part survives and the most general form of v, consistent with the conservation 
laws, is 



X {V a ~ U a ) + {Vj3 - up) 



P. 



(i) 



pRT 



,(!) 



(5) 



where precise form of second order tensor P^l and fictitious velocity u„' need to be determined with the requirements 
of correct i/-Theorem. In subsequent section, I will show that these conditions are fulfilled if 



$ = -\^TTd a logT, 



P {1) - -Jfei T 



dp u a +d a up- —dyUyS a p 



k 2 T d^UySaf 



(6) 



where ki, k 2 and A^ q ' are positive definite scalars related to transport coefficients of the fluid. Here, it needs to be 
noted that the physical picture behind the present model envisions that both pffg and tij, 1 ' are corrections of the 
order of mean free time t . We can relate these new transport coefficients with that of real fluid by performing the 
Chapman-Enskog expansion. The kinetic equation Eq.([3]) along with expression of v given by Eqs.([5]) and are 
the main result of this work. 

It is easy to check that the present formulation gives correct conservation laws for dense gas. This can be seen by 
taking moments of the kinetic equation Eq.Q, which yields time evolution equation for the locally conserved fields 
as 



where 



d t p + d a j a = 0, 

9t jj + d a [p U a Uj + p 5 aJ + (T al ] = 



P 2 ' 

-u +e 



u 2 + e + pj u a + <J a ~< Uy + q a 



(7) 







a =(l + y) cr (K) + P m + P (I) M + ^~ 
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la = <(1 + X) S aP + 
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pRT 
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$+RTu® fdvf 



e- 



D 



(9) 



The kinetic part of the stress tensor and heat flux are conveniently defined in terms of second and third order traceless 
Hermite tensor 



n^ = Ca U 



D + 2 



H 



(2) 

a/3 



1 

D 



(10) 



In terms of these tensors kinetic part of the stress tensor is a^) 



,(K) 



2 R T j dvf ri a a and kinetic part of the heat flux 

is Qa ' = {2RT) AIJ ' J dvf /27ia > ■ We can show that the stress tensor and heat flux has same expression as RET 
if ki, k 2 and are chosen properly. By performing Chapman-Enskog expansion, equation §3§ provides the first 
correction from Maxwell-Boltzmann distribution as 



Cf 



(i) 



2(l + X)WS«a«/J 



V2RT(1 + X )n i a 3) d a logT 



(11) 
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where C is the linearized collision operator. At this stage, we need to provide the exact form of collision operators. 
For example, by using BGK collision operator we obtain 

Ca/3 = -T) (dp u a + d a U(3- ^jdn/UrySapj - k 2 Td 1 u 1 5 al 3, (12) 

where shear viscosity is 77 = (k\ + (1 + x) P) i~ and k 2 r is the bulk viscosity. Similarly for heat flux 

q a = -r^pTX^ +p(l+ X )C p jd a T (13) 

where specific heat at constant pressure C p = (D + 2)/2 R. For BGK model, we have more tuning parameters related 
to transport coefficient then needed. So, without any loss of generality, we can set k\ = 0. So, Prandtl number is 

1-Pr=— ^pTA^Pr (14) 
2pC p (1 + x) 

This means unlike Dufty et al [l5| , in present model we are allowed to set heat conductivity and viscosity coefficients 
independently Thus, if it can be shown that the present model admits //-theorem, it can be claimed that this model 
is first complete phenomenological model for describing dense gas hydrodynamics at the mesoscale level. 

Before proving //-theorem, let us try to understand the physical meaning of Eq. (J5]) . The physical picture behind 
first term on the right hand side of Eq. ([5]) (underlined term A) can be understood in the framework of social force 
model created by Helbing and Molnar in the context of traffic dynamics [l8| . The goal of this term is to keep tagged 
particle separated from other particles, and such crossing of particles can be avoided (only in an average sense) by 
introducing a force which tries to remove them from dense regions. What we want is that at the time r after any 
collision, the probability of having a particle in denser region is smaller. A system mimicking such a motion is 
trajectory of an individual moving on a street trying to avoid being very close to crowded region after time t. Suppose 
he is moving very fast in the region of high density, then he can avoid being in the region of high density at the end of 
time t by moving even more faster. In other words, he experiences a social force which accelerates him in the region of 
high density if he is moving too fast. In the opposite limit, where the individual is moving too slow compared to the 
crowd in dense regions. He can avoid being in the region of high density at the end of time t, if he get decelerated in 
the regions of high density. For hard-spheres such an effect can be modeled by underlined term A in Eq.©. The third 
term C is an hydrodynamic effect and can be understood as microscopic analog of thermophoretic forces [l9| . The 
difference from macroscopic expression is that in stead of transport coefficient a direct dependence on heat velocity 
appears. This may be due to the fact that the current description (Eq. (J3J)) is for short-time motion, whereas transport 
coefficients appear only in long time limit. Similarly, the term B in Eq.© is an hydrodynamic force which reflect the 
tendency to resist locally generated flow field (shear stresses and compression). 

In order to show the existence of //-theorem for the present model, I define the //-function for dense system as 

H (x, t) = J dvf (x, v, t) [log / (x, v, t) - 1] - — ip ^ t] \ (15) 

This choice is similar to that used in proving //-theorem for Enskog equation and is motivated from the work of 
Gremla et al [2(| . Multiplying Eq. §3§ by log / and noting that non- ideal part of entropy is a function of density only, 
we obtain time evolution equation for the H function as 

d t H + J H = J J log fdv + jp ug) 3 a log T + d a u p (16) 



where the flux of //-function (analog of entropy flux) is 

„nid 



Jh = -d a [ '^—u a ) + d a I dv [/(log / - 1) v a ] (17) 



It is interesting to see here that similar to //-function the flux of it also has a contribution totally dependent on 
macroscopic variable p and j. For //-theorem to be valid, entropy production should be positive, which means right 
hand side of Eq. (|16[) must be negative. For Boltzmann collision term or BGK collision term, the first term on the 
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right term is negative. However, this is possible for the last two terms if Eq. ([6]) is valid. Thus, we have proved 
the -ff-theorem. Here, it is interesting to note that the as compared to Boltzmann kinetic theory, the new ingredient 
required to prove if-theorem is just the same as that used in linear irreversible thermodynamics (see for a modern 
perspective [2l|). 

To conclude, in the present manuscript I have presented an alternate mean-field model of hydrodynamics at the 
mesoscale level. Finally, a further interesting simplification of the model happens if we set k\ = &2 = = 0. In 
this limit, using Eq. ([5]) we have 

dt f + v a d a f = - in (P,l T)-f\- d a [J (v a - u a ) x] (18) 
T . ' 

G 

This simplified model can be suitable for numerical implementations in lattice Boltzmann form. In particular, if the 
distribution function / appearing in the term G on the right hand side of Eq. lfTg]) is replaced by / cq , it resembles 
current implementation of hard sphere dynamics in lattice Boltzmann method. Thus, present model provide a self 
consistent framework for hard sphere dynamics in lattice Boltzmann method. Furthermore, similar model can be 
formulated with Boltzmann collision form, which will be very useful for DSMC simulations. 

I am dedicating this work to Dr. I. V. Karlin, teacher, well wisher and friend, in fond memory of the unorthodox 
but comprehensive education in kinetic theory he gave me. I want to thank Professors H. C. Ottinger and Dr. V. 
Kumaran, also my former teachers, for the solid training in non-equilibrium statistical physics and fluid mechanics 
they imparted to me. I am also thankful to Dr. Rochish Thaokar and Dr. S. K. Kwak for the insightful reading of 
this draft. 
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